Non-gaussian statistics from individual pulses of squeezed light 
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We describe the observation of a "degaussification" protocol that maps individual pulses of 
squeezed light onto non-Gaussian states. This effect is obtained by sending a small fraction of 
the squeezed vacuum beam onto an avalanche photodiode, and by conditioning the single-shot ho- 
modyne detection of the remaining state upon the photon-counting events. The experimental data 
provides a clear evidence of phase-dependent non-Gaussian statistics. This protocol is closely related 
to the first step of an entanglement distillation procedure for continuous variables. 
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PACS numbers: 03.67.-a, 42.50.Dv, 03.65.Wj 

Researches on novel schemes to perform quantum key 
distribution (QKD) are presently very active. In that 
field, lots of interest has arisen recently on the use of 
quantum continuous variables (QCV). For instance novel 
QKD schemes using the quadrature components of am- 
plitude and phase modulated coherent states have been 
recently proposed [l| and experimentally demonstrated 
0. It has been shown that such coherent state proto- 
cols are secure against individual gaussian attacks for any 
value of the line transmission [20) an d actually more 
general proofs are presently under study 0, Q . 

An important practical advantage of coherent states 
QKD is that it can in principle reach very high secret bit 
rates Q. However, even in the best possible case, coher- 
ent states QKD will not do much better than photon- 
counting QKD [|| in terms of absolute distance, because 
of the exponential attenuation in optical fibers : at some 
point which is now somewhere between 10 and 100 km, 
one hits a limit where the transmitted secret data gets 
buried into errors of various origins, that range from de- 
tectors dark counts to imperfect data processing. 

In order to qualitatively improve the situation, i.e. to 
go much beyond the attenuation length of a strand of 
fiber, a major challenge is to implement quantum re- 
peaters 0, based upon entanglement distillation and 
(most likely) quantum memories. Ultimately, the se- 
cret qubits would be simply teleported to a remote place, 
with which shared entanglement has been established @ . 
Looking now at entanglement distillation for QCV, a dif- 
ficulty appears quickly : most (if not all) QCV transmis- 
sions so far are using light beams with gaussian statistics. 
However, it has been shown that it is not possible to dis- 
tillate entanglement from agaussian input to a gaussian 
output by gaussian means |9lll0j|. One has to jump "out- 
side" the gaussian domain, though it is possible to reach 
it back at the end, at least in an approximate way pd|. 

In this letter, we experimentally implement a proce- 
dure which we call "degaussification", that maps short 
pulses of squeezed light onto non-Gaussian states. This 
protocol is based upon a post-selection triggered by a 
photon-counting event and uses only simple linear op- 
tical elements. Extending this procedure to entangled 
EPR beams -which is fairly simple in principle- provides 
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FIG. 1: Simplified experimental setup 



the first step of an entanglement distillation procedure 
as proposed in ref. \Ti\ . 

The experimental scheme is presented on Fig. ^ The 
initial pulses are obtained from a titanium-sapphire 
laser (Tiger-CD, Time-Bandwidth Products), delivering 
nearly Fourier-transform limited pulses at 850 nm, with 
a duration of 150 fs, an energy of 40 nJ, and a repeti- 
tion rate of 790 kHz. These pulses are frequency dou- 
bled in a single pass through a thin (100 fiia) crystal of 
potassium niobate (KNbOs), cut and temperature-tuned 
for non-critical type-I phase-matching. The second har- 
monic power is large enough to obtain a significant single- 
pass parametric gain (~ 3 dB) in a similar KNI3O3 crystal 
used in a type-I spatially degenerate configuration. 

Given this relatively high gain, "real" squeezed states 
are actually produced, not only parametric pairs. There- 
fore, higher order terms (beyond pair production) have 
explicitly to be included in the analysis as they play an 
essential role to understand the phase-dependence of the 
data. The detection scheme follows the basic idea of a 
pulsed squeezed light experiment |12|. with two impor- 
tant differences : 

(i) All processing is done in the time domain, not in 
the frequency domain. For each incoming pulse, the bal- 
anced homodyne detection samples one value of the sig- 
nal quadrature in phase with the local oscillator beam 
0- It is then possible to reconstruct the full statistics 
of the signal pulses. The histograms presented below are 
obtained from these individual pulse data. 




FIG. 2: Normalized probability distribution for the (uncon- 
ditioned) squeezed vacuum state, obtained from the pulsed 
homodyne detection. The squeezed quadrature variance is 
1.75 dB below SNL, while the amplified quadrature variance 
is 3.1 dB above. The SNL curve corresponds to the vacuum 
state, where the shot noise variance is taken equal to 1/2. 

(ii) A small fraction (R = 0.115) of the squeezed vac- 
uum beam is taken out from the homodyne detection 
channel. These trigger photons then pass through a spa- 
tial filter (made of two Fourier-conjugated pinholes) and 
a 3 nm spectral filter centered at the laser wavelength, 
before being detected by a silicon avalanche photodiode 
(APD). The detection click is registered simultaneously 
with the homodyne signal, and can be used to post-select 
homodyne events. As we will show, this selection pro- 
vides directly non-gaussian statistics. 

The unconditioned distributions corresponding to the 
squeezed and anti-squeezed quadratures, and to the vac- 
uum noise are plotted on fig. [5] More experimental de- 
tails about the squeezed states generation will be given 
in another publication. The measured squeezing variance 
(with no correction) is 1.75 dB below the shot noise level 
(SNL), in good agreement with the measured deamplifi- 
cation of a probe beam (0.50 or 3 dB) and our evaluation 
of the overall detection efficiency rjtot = f? (1 — R) = 0.66. 
Here 1— R = 0.885 is the transmission of the conditioning 
beamsplitter, and rj = 0.75 is the homodyne detection ef- 
ficiency (see details below). As it can be seen on fig- El the 
experimental data for both quadratures is correctly fit- 
ted by assuming a single-mode parametric gain exp(±2s) 
with s = 0.43, together with the above efficiency r\tot- We 
note however that the deamplification gain of the probe 
beam does not correspond exactly to the inverse of the 
amplification, due to gain-induced-diffraction which dis- 
torts the probe phase fronts |l3j |. Since such multimode 
effects remain reasonably small in our experimental con- 
ditions, we will use the single parameter s to describe 
parametric amplification and deamplification. 

Fig. |3| displays the post-selected output of the homo- 
dyne detection resulting from the degaussification pro- 
tocol, showing a clear dip in the centre of the amplified 
quadrature distribution. The theoretical curves repre- 
sented on the same figure are obtained from a simple 
single-mode model detailed below. This model takes into 
account the measured parametric gain, together with var- 
ious experimental imperfections (losses, imperfect mode- 
matching, electronic noise, dark counts and modal purity, 
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FIG. 3: Experimental (dots) and theoretical (line) quadrature 
distribution of the post-selected homodyne measurements for 
the amplified quadrature (a) and the squeezed one (b), nor- 
malized as in fig. Parameters used in the calculation are 
s = 0.43, R = 0.115, v = 0-75 and £ = 0.7 . 

see below for details) , and it is clearly in good agreement 
with the experimental data. 

The origin of the observed effect can be analyzed in 
different ways. A first insight can be obtained by con- 
sidering the homodyne detection of a conditional single 
photon state, observed in 01 . In this experiment, the 
authors separate the two photons from a parametric pair, 
and one of them is used as a trigger on a photon counter, 
while the other one is sent to an homodyne detection. 
In the ideal case, such an experiment would measure the 
probability density P(x) of the n = 1 Fock state, which 
is non-gaussian since P(0) = 0. Though this experiment 
provides a first idea of the origin of the non-gaussian fea- 
tures, it is not enough to explain our observations. Actu- 
ally, we see phase-dependent effects (while a n = 1 Fock 
state is phase-independent), and in our set-up there is no 
explicit separation of the photon pair. 

We have carried out a calculation taking into account 
the expansion of the squeezed state in a Fock state basis, 
including terms up to n = 10, which is enough for our 
degree of squeezing. The calculation is done for an arbi- 
trary value of the conditioning beamsplitter reflectivity, 
and takes into account the various imperfections of the 
experiment. This calculation is straightforward but te- 
dious, and one can actually get a good physical insight 
by considering the restricted simple case of an expansion 
of the squeezed vacuum up to n — 4, and a beamsplitter 
reflectivity r = -c 1. The squeezed vacuum | W s ) 
can then be written as : 

|* a > = a|0)+/3|2>+ 7 |4) (1) 

With our degree of squeezin g s = 0.43, one has a = 
0.96, = 0.27 and 7 = 0.10 [L^. This state then gets 
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FIG. 4: Phase-dependent quadrature distributions of the con- 
ditioned homodyne measurements, together with the vacuum 
reference (line and dots) . The thick solid line is obtained from 
eq. with R — 0.01. The thin gray line is obtained from 
the complete calculation and R — 0.115. Fig. (a) corresponds 
to the amplified quadrature while fig. (b) shows the squeezed 
one. The squeezing parameter is s — 0.43, and perfect single 
mode detection efficiency has been assumed. 

mixed with the vacuum at the beamsplitter, resulting 
in a two-mode entangled squeezed state. Denoting as 
r, t the reflectivity and transmittance of the beamsplitter 
(r 2 + 1 2 = 1), the output state is : 

I 9 a , mt ) = (a I 0)i +i 2 /? | 2)i + i 4 7 I 4)i) | 0) 2 

+ (V2rt/3 | l)i + 2rt 3 7 | 3)i) | l)a + 0(2) (2) 

where | .)i denotes the state sent to the homodyne de- 
tection, while | .) 2 stands for the state sent to the APD. 
The term 0(2) denotes Fock state terms higher than 1 on 
the APD beam, which will be neglected in this simplified 
calculation, given our assumption r <C 1. Finally, post- 
triggering on the APD photon-counting events reduces 
the state detected by the homodyne detection to : 

I * co , id ) « | 1) + V2 7 i 2 | 3) (3) 

The prediction of this calculation is shown on fig. 0] As 
it could be expected, we do obtain phase-dependent non- 
gaussian statistics. These features are related to high 
order terms beyond pair production which play an essen- 
tial role in our analysis. 

In this simplified calculation we have assumed r << 1, 
and the predicted dip in the centre of the probability 
distribution goes down to zero. When the beamsplitter 
reflectivity is increased, Fock state terms with n > 1 may 
no longer be neglected on the APD beam, and the central 
dip has a non-zero value. Strictly speaking, this is not an 
experimental imperfection, but an intrinsic feature of the 
conditioned state for larger R, which clearly appears on 
the result of the full calculation also displayed on Fig. 4. 

In order to characterize experimental imperfections, 
let us emphasize that the homodyne detection and the 
photon-counting detection have quite different draw- 
backs. The homodyne detection is not sensitive to "real" 
photons that are in modes unmatched with the detected 
(local oscillator) mode, but it is quite sensitive to vac- 
uum modes which couple into this detected mode. On 
the other hand, the photon-counting detection is not sen- 
sitive to vacuum noise, but it will detect photons in any 



modes. Correspondingly, two experimental parameters 
must be used : an homodyne efficiency parameter rj, 
which measures the overlap between the desired signal 
mode and the detected mode |l5j ; and a modal pu- 
rity parameter £, which characterizes which fraction of 
the detected photons are actually in the desired signal 
mode |l6j . In the simplest approach, the homodyne effi- 
ciency can be modelized by a lossy beamsplitter, taking 
out desired correlated photons. On the other hand, the 
modal purity £ in our experiment cannot be modelized 
by another lossy beamsplitter, because a small value of 
£ corresponds to unwanted firings of the APD, for which 
a squeezed vacuum is still measured at the homodyne 
detection port. More precisely, the measured proba- 
bility distribution for a quadrature x will be taken as 

P(x) = £ P CO nd(x) + (1 - £) Puncond(x), where P C ond(x) 

and P U ncond{x) are respectively the conditioned and un- 
conditioned probability distributions, which depend on 
the values of s, R and rj. 

It is then easy to determine values of the parameters 
rj and £ fitting the experimental data. The procedure 
to measure rj is well established from squeezing exper- 
iments [l2j, and it can be cross-checked by comparing 
the classical parametric gain and the measured degree of 
squeezing. The procedure to measure £ is less usual, and 
amounts to evaluate how many unwanted photons make 
their way through the spatial and spectral filters which 
are used on the photon counting channel. Ultimately, 
this estimated value of £ must fit with the observed con- 
ditional probability distribution, since r/ is independantly 
obtained from squeezing measurements. 

Experimentally, this procedure turns out to be quite 
successful, and for instance we have plotted on fig. |3| 
the amplified and deamplified conditional probability 
distributions, using as parameters the parametric gain 
exp(2s) = 2.36, the homodyne efficiency rj = 0.75, and 
the modal purity parameter £ = 0.7. We note that the 
value of s is evaluated from the measured squeezing (see 
fig. |2J), while n is obtained as rj = rjTi]j { i]D, where the 
overall transmission ryr = 0.94, the mode-matching visi- 
bility r]H = 0.92, and the detectors efficiency r\n = 0.945 
are independantly measured. Finally, the modal purity 
£ is fitted to the data, and cross-checked as the ratio 
between the expected and actual APD counting rates. 

In a last step, we have analysed our data using the 
standard techniques of quantum tomography. We have 
recorded an histogram with 40 bins for 6 different quadra- 
ture phase values 9, and about 5000 points for each his- 
togram were acquired in a 3 hours experimental run. 
The Wigner function displayed on fig. [5] was then recon- 
structed using the Radon transform \V?\. applied to the 
symetrized experimental data (P(xg)+P(— xg))/2, with- 
out any correction for measurement efficiency. It shows a 
clear dip at the origin, with a central value of 0.067 while 
the maximum is at 0.12. 

As usual, the conditions to get negative values of the 
measured Wigner function are rather stringent, and re- 
quire the presence of a dip into the distribution proba- 




FIG. 5: (a) Theoretical Wigner function W of the output 
state of the "degaussification" protocol, assuming s = 0.43, 
R = 0.115 and perfect detection (rj — £ = 1). (b) Recon- 
structed Wigner function from the experimental data (77 = 
0.75, £ = 0.7). The values of W at the origin of phase space 
are respectively W th (0,0) = -0.26, and W eX p{0,0) = 0.067. 



bility associated to the squeezed quadrature. Given our 
experimental parameters, this requires a modal purity £ 
better than 0.85, which was not experimentally attain- 
able while keeping the APD count rate above a few tens 
per second. Nevertheless, we point out that by correct- 
ing for the homodyne efficiency, the evaluated Wigner 
function of the prepared state (just before homodyne 
detection) does assume a negative value at the origin, 
W cor (0,0) « —0.06. Another interesting feature is that 
the non-gaussian dip on the amplified quadrature is quite 
robust to losses and therefore can be easily observed with 
our experimental parameters. This is associated with a 
similarly robust "squeezed volcano shape" of the Wigner 
function. 

We have described the first experimental observation 
of a "degaussification" protocol, mapping individual fem- 
tosecond pulses of squeezed light onto non-Gaussian 
states, by using only linear optical elements and an 
avalanche photodiode. The observed effect is closely re- 
lated to the first step of an entanglement distillation pro- 
cedure for gaussian quantum continuous variables 
This work should contribute to the future development 
of quantum repeaters and long-range quantum cryptog- 
raphy using continuous variables entanglement. 
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